Uncertainties in nuclear transition matrix elements for neutrinoless f3/3 decay II: the 

heavy Majorana neutrino mass mechanism 
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Employing four different parametrizations of the pairing plus multipolar type of effective two- 
body interaction and three different parametrizations of Jastrow-type of short range correlations, 
the uncertainties in the nuclear transition matrix elements AiS"' due to the exchange of heavy 



Majorana neutrino for the + — > + transition of neutrinoless double beta decay of 94 Zr, 96 Zr, 
98 Mo, 100 Mo, 
be around 35% 

parametrization of Jastrow short range correlations, the uncertainties are found to be smaller than 
20%. 



104 Ru, 110 Pd, 128 . 130 Te and 150 Nd isotopes in the PHFB model are estimated to 
Excluding the nuclear transition matrix elements calculated with Miller-Spenser 



PACS numbers: 21.60.-n, 23.40.-s, 23.40.Hc 



o 

=5 



(N 
> 
O 

in 
o 



I. INTRODUCTION 

In addition to establishing the Dirac or Majorana na- 
ture of neutrinos, the observation of (Pf3) 0l/ decay is a 
convenient tool to test the lepton number conservation, 
possible hierarchies in the neutrino mass spectrum, the 
origin of neutrino mass and CP violation in the leptonic 
sector. Further, it can also ascertain the role of various 
gauge models associated with all possible mechanisms, 
namely the exchange of light neutrinos, heavy neutri- 
nos, the right handed currents in the left-right symmet- 
ric model (LRSM), the exchange of sleptons, neutrali- 
nos, squarks and gluinos in the R p -violating minimal 
super symmetric standard model, the exchange of lep- 
toquarks, existence of heavy sterile neutrinos, compos- 
iteness, extradimensional scenarios and Majoron models, 
allowing the occurrence of (/3/3) 0j/ decay. Stringent limits 
on the associated parameters have already been extracted 
from the observed experimental limits on the half-life of 
(P~l3~) 0l/ decay [l| and presently, all the experimental 
attempts are directed for its observation. The experi- 
mental and theoretical studies devoted to (f3f3) 0v decay 
over the past decades have been recently reviewed by 
Avignone et al. |2| and references there in. 

Presently, there is an increased interest to calculate re- 
liable NTMEs for (/3~ f3~ ) Qv decay due to the exchange 
of heavy Majorana neutrinos, in order to ascertain the 
dominant mechanism contributing to it The lep- 

ton number violating [f}~ ) decay has been studied 
by Vergados by taking a Lagrangian consisting of left- 
handed as well as right-handed leptonic currents [5|. In 
the QRPA, the (0~0~) Ov decay due to the exchange of 
heavy Majorana neutrinos has been studied by Tomoda 
@. The decay rate of (/3~ ' P~) Qv mode in the LRSM has 
been derived by Doi and Kotani Q- Hirsch et al. [1] have 
calculated all the required nuclear transition matrix ele- 



ments (NTMEs) in the QRPA and limits on the effective 
light neutrino mass heavy neutrino mass (Mat), 

right handed heavy neutrino (Mr), (A), (rj) and mix- 
ing angle tan£ have been obtained. The heavy neutrino 
mechanism has also been studied in the QRPA without 
and with pn- pairing jioj . In the heavy Majorana neu- 
trino mass mechanism, Simkovic et al. have studied 
the role of induced weak magnetism and pseudoscalar 
terms and it was found that they are quite important 
in 48 Ca nucleus. The importance of the same induced 
currents in both light and heavy Majorana neutrino ex- 
change mechanism has also been studied using the pn- 
RQRPA [3 as well as SRQRPA 0. 

In spite of the remarkable success of the large scale 
shell model (LSSM) calculations of Strassbourg-Madrid 
group [l3| , there is a necessity of large configuration mix- 
ing to reproduce the structural complexity of medium 
and heavy mass nuclei. On the other hand, the QRPA 
and its extensions have emerged as successful models by 
including a large number of basis states and in correlat- 
ing the single-/? GT strengths and half-lives of (/3~/3 _ )2i/ 
decay in addition to explaining the observed suppression 
of M 2 „ HIE!. In the mass region 90 < A < 150, there 
is a subtle interplay of pairing and quadrupolar correla- 
tions and their effects on the NTMEs of (/3^/3~) 0v decay 
have been studied in the interac ting shell model (ISM) 
mm> deformed QRPA model andprojected- 
Hartree-Fock-Bogoliubov (PHFB) model [22I. [23l| . 

The possibility to constrain the values of the gauge 
parameters using the measured lower limits on the 
(/3 _ /3~) „ decay half-lives relies heavily on the model de- 
pendent NTMEs. Different predictions are obtained by 
employing different nuclear models, and within a given 
model, varying the model space, single particle energies 
(SPEs) and effective two-body interaction. In addition, 
a number of issues regarding the structure of NTMEs, 
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namely the effect of pseudoscalar and weak magnetism 
terms on the Fermi, Gamow- Teller and tensorial NTMEs 
[13, H5[, the role of finite size of nucleons (FNS) as well as 
short range correlations (SRC) vis-a-vis the radial evo- 
lution of NTMEs [H, f26M28| and the value of the axial- 
vector coupling constant gA are also the sources of un- 
certainties and remain to be investigated. 

It was observed by Vogel [2!| that in case of well stud- 
ied 76 Ge, the calculated decay rates T^ 2 differ by a factor 
of 6-7 and consequently, the uncertainty in the effective 
neutrino mass (w„) is about 2 to 3. Thus, the spread be- 
tween the calculated NTMEs can be used as the measure 
of the theoretical uncertainty. In case the (/3/3) „ deca y of 
different nuclei will be observed, Bilenky and Grifols [30] 
have suggested that the results of calculations of NTMEs 
of the (/3~/3~) 0v decay can be checked by comparing the 
calculated ratios of the corresponding NTMEs-squared 
with the experimentally observed values. 

Bahcall et al. [31] and Avignone et al. [32[ have cal- 
culated averages of all the available NTMEs, and their 
standard deviation is taken as the measure of theoretical 
uncertainty. On the other hand, Rodin et al. (33| have 
calculated nine NTMEs with three sets of basis states and 
three realistic two-body effective interactions of charge 
dependent Bonn, Argonne and Nijmen potentials in the 
QRPA as well as RQRPA and estimated the theoretical 
uncertainties by making a statistical analysis. It was no- 
ticed that the variances are substantially smaller than the 
average values and the results of QRPA, albeit slightly 
larger, are quite close to the RQRPA values. Faessler and 
coworkers have further studied uncertainties in NTMEs 
due to short range correlations using unitary correlation 
operator method (UCOM) [2H and self-consistent cou- 
pled cluster method (CCM) [27[. 

The PHFB model has the advantage of treating the 
pairing and deformation degrees of freedom on equal foot- 
ing and projecting out states with good angular momen- 
tum. However, the single j3 decay rates and the distri- 
bution of GT strength, which require the structure of 
the intermediate odd Z-odd N nuclei, can not be stud- 
ied in the present version of the PHFB model. In spite 
of this limitation, the PHFB model in conjunction with 
pairing plus quadrupole-quadrupole (PQQ) [34[ has been 
successfully applied to reproduce the lowest yrast states, 
electromagnetic properties of the parent and d aug hter 
nuclei, 
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and the measured (/3~ 0~) 2v decay rates 
In the PHFB formalism, the existence of an inverse corre- 
lation between the quadrupole deformation and the size 
of NTMEs M 2v , M (0l/) and M§ v) has been observed 
[H HI]. Further, it has been noticed that the NTMEs 
are usually large for a pair of spherical nuclei, almost con- 
stant for small deformation, suppressed depending on the 
difference in the deformation A/?2 of parent and daughter 
nuclei and having a well defined maximum when A/?2 = 

mm. 

In Ref. 



37j , a statistical analysis was performed for 



neutrino, calculated in the PHFB model with four dif- 
ferent parameterizations of pairing plus multipolar type 
of effective two-body interaction [23| and two (three) dif- 
ferent parametrization of Jastrow type of SRC [271 ] . In 
confirmation with the observation made by Simkovic et 
al. [13]' ^ was noticed that the Miller-Spenser type of 
parametrization is a major source of uncertainty and 
its exclusion reduces the uncertainties from 10%-15% 
to 4%-14%. Presently, the same procedure has been 
adopted to estimate the theoretical uncertainties asso- 
ciated with the NTMEs M§ v) for {P~P') Qv decay due 
to the exchange of heavy Majorana neutrino. In Sec. 
II, a brief discussion of the theoretical formalism is pre- 
sented. The results for different parameterizations of the 
two-body interaction and SRC vis-a-vis radial evolution 
of NTMEs are discussed in Sec III. In the same section, 
the averages as well as standard deviations are calculated 
for estimating the theoretical uncertainties. Finally, the 
conclusions are given in Sec. IV. 



II. THEORETICAL FORMALISM 

In the charged current weak processes, the current- 
current interaction under the assumption of zero mass 
neutrinos leads to terms which, except for vector and 
axial vector parts, are proportional to the lepton mass 
squared, and hence negli gibl e. However, it has been re- 
ported by Simkovic et al. |24Ll25| that the contribution of 
the pseudoscalar term is equivalent to a modification of 
the axial vector current due to PC AC and greater than 
the vector current. The contributions of pseudoscalar 
and weak magnetism terms in the mass mechanism can 
change M^ ^ upto 30% and the chan ge i n M^ v ^ is con- 
siderably larger. In the shell- model [lj, E0J, IBM [H 
and GCM+PNAMP [|oj|, the contributions of these pseu- 
doscalar and weak magnetism terms to have been 
also investigated. However, it has been reported by Suho- 
nen and Civitarese [4l| that these contributions are rel- 
atively small and can be safely neglected. Therefore, the 
investigation of this issue is of definite interest and is 
reported in the present work. 

In the two nuclcon mechanism, the half-life T®J 2 for 

the + — > + transition of {/3~/3~) 0u decay due to the 
exchange of heavy Majorana neutrino between nucleons 
having finite size is given by @, 0] 
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extracting uncertainties in eight (twelve) NTMEs for 
(fl~j3~) 0l , decay due to the exchange of light Majorana 



where m p is the proton mass and 

(Mn)- 1 = V U 2 el mr\ m, > 1 GeV, (2) 

z — '1 

and in the closure approximation, the NTMEs M^'^ is 
of the form [H M, Wk 
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The exchange of heavy Majorana neutrinos gives rise 
to short ranged neutrino potentials, which with the con- 
sideration of FNS are given by 



where 
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where f a h (qr n m) = jo {qrnm) for a = F as well as GT 
and f T h {qr inn ) = h (q r urn J • 

Further, the h F (q), hQp(q) and hr{q) are written as 
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where the form factors are given by 
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with gv — 1.0, 3a = 1-254, k = [i p — (x n = 3.70, Ay = 
0.850 GeV, A^ = 1.086 GeV and m w is the pion mass. 

Substituting Eq. ©-Eq. (TJT]) in Eq. ©, there is 
one term, associated with hp, Eq. ([§]), contributing to 
Mph, while Mark has four terms, denoted by Mgt-aa, 
Mqt_ ap, Mqt-pp and Mqp_mm-, which correspond to 
the four terms in her, Eq. (|10j) . The tensor contribution, 
AfTh, has three terms, denoted by Mp-AP, Mp~pp and 
Mp-MM, which correspond to the three terms in hp, 
Eq. (fTTj) . Their contributions to the total nuclear matrix 
element are discussed in Sec. III. 

The short range correlations (SRC) arise mainly from 
the repulsive nucleon-nucleon potential due to the ex- 
change of p and u> mesons and have been incorporated by 



using effective transition operator [42j | , the exchange of oj- 
meson @, UCOM HQ and the self-consistent CCM 
[27] |. The SRC can also be incorporated phenomenologi- 
cally by Jastrow typ e of correlations with Miller-Spenser 
parametrization [45J. Further, it has been shown in the 
self-consistent CMM [27[ that the SRC effects of Argonne 
and CD-Bonn two nucleon potentials are weak and it is 
possible to parametrize them by Jastrow type of correla- 
tions within a few percent accuracy. Explicitly, 
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where a = 1.1, 1.59 and 1.52 fm- 2 , b = 0.68, 1.45 and 
1.88 fm- 2 and c = 1.0, 0.92 and 0.46 for Miller-Spencer 
parametrization, CD-Bonn and Argonne V18 NN poten- 
tials, respectively. In this work the NTMEs Afjy are 
calculated in the PHFB model for the above mentioned 
three sets of parameters for the SRC, denoted as SRC1, 
SRC2 and SRC3, respectively. 

In Fig.l, we plot the neutrino potential Hn (r, A)= 
Hp h {r, A) /(r) with the three different parametrizations 
of SRC. It is noticed, that the potentials due to FNS and 
FNS+SRC3 are peaked at the origin where as the peaks 
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FIG. 1: Radial dependence of H N (r, A)= H Fh (r, A) f(r) for 
the three different parameterizations of the SRC. In the case 
of FNS, f{r) = 1. 



due to FNS+SRC1 and FNS+SRC3 are at r w 0.6 fm and 
r 0.5 fm, respectively. The shapes of these functions 
have definite influence on the radial evolution of NTMEs 
M^ v ^ for (/3~/3~) 0l/ decay due to the exchange of heavy 
Majorana neutrino as discussed in Sec. III. 



The calculation of M^ u ^ in the PHFB model has been 
discussed in our earlier work [22], |37j and one obtains 
the following expression for NTMEs Mf u) of (/?" /3") 0i/ 
decay [37| 
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and the expressions for calculating rc. , 
n (z,7V)Xz.+2jy-2)(^), /z,jv and F z ,n(0) are given in 

Refs. [II Eg. 

The calculation of matrices fz.N an d Fz,n{9) re- 
quires the amplitudes (ui m ,Vi m ) and expansion coeffi- 
cients Cij tm , which specify the axially symmetric HFB 
intrinsic state |$o) with K = 0. Presently, they are 
obtained by carrying out the HFB calculations through 
the minimization of the expectation value of the effective 



Hamiltonian given by (23j 

H = H sp + V(P) + V(QQ) + V{HH) (15) 

where H Bp , V(P), V(QQ) and V(HH) denote the 
single particle Hamiltonian, the pairing, quadrupole- 
quadrupole and hexadecapole-hexadecapole part of the 
effective two-body interaction, respectively. The HH 
part of the effective interaction V(HH) is written as [23[ 



V(HH) = - (f ) E £(-l)>|r%„(M)l7>0Vn,-v(MM o+oj a 5 a, (16) 

i 



with X4 = 0.2442 X2^4~ 2/3 &~ 4 for T = 1, and twice of 
this value for T = case, following Bohr and Mottelson 

In Refs. [H, [H, the strengths of the like par- 
ticle components Xvv anc ^ Xnn of the QQ interaction 
were kept fixed. The strength of proton-neutron (pn) 
component x-pn was varied so as to reproduce the ex- 



citation energy of the 2 + state E 2 + for the consid- 
ered nuclei, namely 94 ' 96 Zr, 94,96,98,ioo MOj 98 ,ioo,io4 RUj 

104,U0 Pd| 110 Cd; 128,130^ 128,130 Xej 150 Nd and 150 Sm 

as closely possible to the experimental values. This is de- 
noted as PQQl parametrization. Alternatively, one can 
employ a different parametrization of the X2pn, namely 
PQQ2 by taking X2 PP = X2nn = X2 pn /2 and the ex- 
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citation energy E 2 + can be reproduced by varying the 
X2pp- Adding the HH part of the two-body interaction to 
PQQl and PQQ2 and by repeating the calculations, two 
more parameterizations of the effective two-body interac- 
tions, namely PQQHHl and PQQHH2 were obtained 

The four different parameterizations of the effective 
pairing plus multipolar correlations provide us four dif- 
ferent sets of wave functions. With three different param- 
eterizations of Jastrow type of SRC and four sets of wave 
functions, sets of twelve NTMEs Afjy are calculated 
for estimating the associated uncertainties in the present 
work. The uncertainties associated with the NTMEs 
M^' for (/3 - /3 - ) „ decay are estimated statistically by 
calculating the mean and the standard deviation defined 
by 
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III. RESULTS AND DISCUSSIONS 



The model space, SPE's, parameters of PQQ type of 
effective two-body interactions and the method to fix 
them have already been given in Refs. [22|, [3|1 |36| . 
It turns out that with PQQl and PQQ2 parameteri- 
zations, the experimental excitation energies of the 2 + 
state E 2 + [47]] can be reproduced within about 2% ac- 
curacy. The electromagnetic properties, namely reduced 
B(E2:0 + —> 2 + ) transition probabilities, deformation pa- 
rameters /32, static quadrupole moments Q(2 + ) and gy- 
romagnetic factors g(2~^) are in overall agreement with 
the experimental data [HI, |4^] . 



Short range correlations and radial evolutions 
of NTMEs 



TABLE I: Calculated NTMEs M# v) in the PHFB model with 
four different parameterization of effective two-body interac- 
tion and three different parameterizations of Jastrow type of 
SRC for the (/T/T)^ decay of 94 ' 96 Zr, 98 > 100 Mo, 104 Ru, 
110 Pd, 128 ' 130 Te and 150 Nd isotopes due to the exchange of 
heavy Majorana neutrino exchange, (a), (b), (c) and (d) de- 
note PQQl, PQQHHl, PQQ2 and PQQHH2 parameteri- 
zations, respectively. See the footnote in page 3 of Ref. [33] 
for further details. 



Nuclei 




F 




F+S 










SRC1 


SRC2 


SRC3 


y4 Zr 


(a.) 


236.9498 


77.5817 


138.2606 


191.3897 




(b) 


220.3794 


72.4285 


128.7496 


178.0783 




(c) 


205.8370 


72.9303 


124.3248 


168.5705 




(d) 


211.0437 


68.9323 


122.9710 


170.3572 


96 Zr 


(a) 


177.7479 


56.4909 


102.4434 


142.8831 




(b) 


185.5251 


59.5338 


107.2877 


149.3117 




(c) 


170.8199 


54.2382 


98.4051 


137.2870 




(d) 


175.4730 


56.0746 


101.2963 


141.1240 


98 Mo 


(a) 


355.1915 


117.0804 


208.2494 


287.5615 




(b) 


346.1118 


116.4967 


204.5667 


281.0515 




(c) 


358.5109 


118.0563 


210.1150 


290.2080 




(d) 


343.4160 


115.2077 


202.6977 


278.7158 


100 Mo 


(a) 


365.8004 


122.2000 


215.8882 


296.9869 




(b) 


361.9877 


122.6611 


214.7455 


294.4297 




(c) 


368.4056 


123.2364 


217.5391 


299.1598 




(d) 


328.9795 


111.4464 


195.1601 


267.5869 


104 Ru 


(a) 


274.0700 


89.7666 


160.7925 


222.1151 




(b) 


264.9015 


88.1515 


156.2893 


215.1076 




(c) 


258.2796 


84.6746 


151.6002 


209.3600 




(d) 


247.0603 


82.3208 


145.8435 


200.6645 


uo pd 


(a) 


424.6601 


140.3359 


249.6835 


344.3187 




(b) 


379.9404 


127.4915 


224.6563 


308.6907 




(c) 


407.2163 


134.6824 


239.4733 


330.1888 




(d) 


390.3539 


130.6314 


230.5392 


316.9996 


128 Te 


(a) 


190.5325 


62.4373 


111.5143 


154.1796 




(b) 


231.8024 


77.4559 


136.7936 


188.1893 




(c) 


220.7156 


73.5158 


130.0810 


179.0960 




(d) 


235.4814 


78.6367 


138.9052 


191.1366 


130 Te 


(a) 


236.0701 


81.5493 


141.3447 


192.7610 




(b) 


231.5921 


79.3844 


138.1901 


188.8492 




(c) 


233.0024 


80.4020 


139.4400 


190.2194 




(d) 


230.5282 


78.9888 


137.5307 


187.9675 


150 Nd 


(a) 


163.8037 


55.8968 


97.8169 


133.6912 




(b) 


130.1364 


43.8840 


77.3178 


105.9993 




(c) 


160.2720 


54.6713 


95.6942 


130.8005 




(d) 


131.9781 


44.6741 


78.5433 


107.5715 



In the approximation of finite size of nucleons with 
dipole form factor (F) and finite size plus SRC (F+S), 
the theoretically calculated twelve NTMEs M^ v ' us- 
ing the four sets of HFB wave functions generated with 
PQQl, PQQHHl, PQQ2 and PQQHH2 parameteri- 
zations of the effective two-body interaction and three 
different parameterizations of Jastrow type of SRC for 
94 < 96 Zr, 98 < 100 Mo, 104 Ru, 110 Pd, i28d30 Te and i50 Nd igo _ 

topes are presented in Table Q] 

To analyze the role of different components of NTME 
M^ v \ the decomposition of the latter into Fermi, differ- 
ent terms of Gamow- Teller and tensor matrix elements of 



100 Mo are presented in Table [TTI for PQQl parametriza- 
tion. From the inspection of Table HH the following ob- 
servations emerge. 

(i) The contribution of conventional Fermi matrix el- 
ements Mph = Mp-vv is about 20% to the total 
matrix element. 

(ii) The Gamow- Teller matrix element is noticeably 
modified by the inclusion of the pseudoscalar and 
weak magnetism terms in the hadronic currents. 
While Mct-pp increases the absolute value of 
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TABLE II: Decomposition of NTMEs for the {P~P~) 0v decay 
of 100 Mo including finite size effect (F) and SRC (F+S) for 
the PQQ1 parametrization. 



NTMF.s 


p 








SRC1 


SRC2 


SRC3 




UO.UiiO 






(\A OKI c 


Mgt-aa 


-370.5960 


-144.5650 


-242.3340 


-316.3250 


Mqt-ap 


174.4640 


43.3631 


93.9737 


137.5100 


Mgt-pp 


-66.3082 


-8.3767 


-28.5936 


-48.0727 


Mgt-mai 


-41.7693 


16.3949 


7.6213 


-13.3421 


Mgt 


-304.2095 


-93.1837 


-169.3326 


-240.2298 


Mt-ap 


9.4369 


9.2393 


10.0332 


10.0610 


Mt—pp 


-3.6622 


-3.5528 


-3.9226 


-3.9386 


Mt-mm 


1.2567 


1.1163 


1.3438 


1.3722 


M t 


7.0314 


6.8028 


7.4544 


7.4945 



Ml '\ 365.8004 122.2000 215.8882 296.9869 



600.0 n 




12 3 4 

r (fm) 

FIG. 2: Radial dependence of C^" ] (r) for the (P~P~) 0v de- 
cay of 100 Mo isotope. 



Mqt-aa, Mqt-ap has a significant contribution 
with opposite sign in all cases. The term Mqt-mm 
is smaller than others, and the introduction of short 
range correlations changes its sign. 

(iii) The tensor matrix elements have a very small con- 
tribution, smaller than 2%, to the total transition 
matrix elements. 

(iv) The inclusion of short range correlations changes 
the nuclear matrix elements significantly, whose ef- 
fects are large for the Gamow- Teller and Fermi ma- 
trix elements but small in the case of tensor ones. 

(v) The Miller- Spencer parameterization of the short 
range correlations, SRC1, cancels out a large part 
of the radial function H^, as shown in Fig. 1. The 
same cancellation reduces the calculated matrix el- 
ements to about one third of its original value. The 
other two parameterizations of the short range cor- 
relations, namely SRC2 and SRC3, have a sizable 
effect, which is in all cases much smaller than SRC1. 

With respect to the point nucleon case, the change 
in M^ v) is about 30%-34% due to the FNS. With the 
inclusion of effects due to FNS and SRC, the NTMEs 
change by about 75%-79%, 58%-62% and 43%-47% for 
F+SRC1, F+SRC2 and F+SRC3, respectively. It is 
noteworthy that the SRC3 has practically negligible ef- 
fect on the finite size case. Further, the maximum varia- 
tion in M§ v) due to PQQHH1, PQQ2 and PQQHH2 
parametrization with respect to PQQ1 interaction are 
about 24%, 18% and 26% respectively. 

In the QRPA [H,[27j], ISM [1|| and PHFB [MUzI, the 
radial evolution of M^ 0v ' due to the exchange of light 
Majorana neutrino has already been studied. In both 
QRPA and ISM calculations, it has been established that 
the contributions of decaying pairs coupled to J = and 



J > almost cancel beyond r «3 fm and the magnitude 
of C' 0l/ ) for all nuclei undergoing (/3 _ /3~) „ decay have 
their maximum at about the internucleon distance r ~1 
fm. These observations were also made in the PHFB 
model [28|, |37|. Similarly, the radial evolution of Afjy 
can be studied by defining 

mT ] = JcT } (r)dr (19) 

The radial evolution of M^"' has been studied for four 
cases, namely F, F+SRC1, F+SRC2 and F+SRC3. To 
make the effects of finite size and SRC more transparent, 
we plot them for 100 Mo in Fig. 2. In case of finite sized 
nucleons, the Cjy are peaked at r «0.5 fm and with the 
addition of SRCl and SRC2, the peak shifts to about 0.8 
fm. However, the position of peak is shifted to 0.7 fm 
for SRC3. In Fig. 3, we plot the radial dependence of 
C [ ° u) for six nuclei, namely 96 Zr, 100 Mo, 110 Pd, 128 > 130 Te 
and 150 Nd and the same observations remain valid. Also, 
the same features in the radial distribution of cjy 1 ^ are 
noticed in the cases of PQQ2, PQQHH1 and PQQHH2 
parametrizations. 



B. Uncertainties in NTMEs 

The uncertainties associated with the NTMEs M^ u ^ 
for {j3~/3~) 0u decay are estimated by preforming a sta- 
tistical analysis by using Eas. fTTl) and (TT51) . In Table |H 

sets of twelve NTMEs M { ° u) of 94 - 96 Zr, 98 - 100 Mo, 110 Pd, 
128 ' 130 Te and 150 Nd isotopes are displayed, which are em- 
ployed to calculate the average values as well as 

uncertainties AM^ tabulated in Table Hill for the bare 
axial vector coupling constant gA = 1.254 and quenched 
value of gA — 1-0. 
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FIG. 3: Radial dependence of C$ v) (r) for the {P~P~) 0v decay of 96 Zr, 100 Mo, 110 Pd, 12; 
Fig., (a), (b), (c) and (d) correspond to F, F+SRC1, F+SRC2 and F+SRC3, respectively. 



3 Te and 150 Nd isotopes. In this 



It turns out that in all cases, the uncertainties AM 
are about 35% for gA = 1.254 and gA = 1-0. Further, we 

estimate the uncertainties for eight NTMEs calcu- 
lated using the SRC2, and SRC3 parameterizations and 
the uncertainties in NTMEs reduce to about 16% to 20% 
with the exclusion of Miller-Spenser type of parametriza- 
tion. In Table IIV| average NTMEs for case II along with 
NTMEs calculated in other models have been presented. 
It is noteworthy that in the models employed in Refs. 
0, IS IH, effects due to higher order currents have not 
been included. We also extract lower limits on the ef- 
fective mass of heavy Majorana neutrino (Mjv) from the 
largest observed limits on half-lives T®" 2 of (j3~ j3~) 0v de- 
cay. The extracted limits are (M N ) > 5.67lo'j 4 x 10 7 
GeV and > 4.06±£^ xlO 7 GeV, from the limit on half- 
life T°J 2 >3.0xl0 24 yr of 130 Te [H for g A = 1.254 and 
gA — 1.0, respectively. 



IV. CONCLUSIONS 

We have employed the PHFB model, with four differ- 
ent parameterizations of pairing plus multipole effective 
two body interaction, to generate sets of four HFB in- 
trinsic wave functions, which reasonably reproduced the 
observed spectroscopic properties, namely the yrast spec- 
tra, reduced B(E2:0 + — » 2 + ) transition probabilities, 
static quadrupole moments Q(2 + ) and g-factors <?(2 + ) 
of participating nuclei in (/?~ /3~) 2 „ decay, as well as 
their Mi v (35l . l36l | . Considering three different parame- 
terizations of Jastrow type of SRC, sets of twelve NTMEs 
M { ° u) for the study (/3~ 0~) o decay of 94 - 96 Zr, 98 > 100 Mo, 
104 Ru, 110 Pd, 128 > 130 Te and 150 Nd isotopes in the heavy 
Majorana neutrino mass mechanism have been calcu- 
lated. 

The study of effects due to finite size of nucleons and 
SRC reveal that in the case of heavy Majorana neutrino 
exchange, the NTMEs change by about 30%-34% due to 
finite size of nucleons and the SRC1, SRC2 and SRC3 
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TABLE III: Average NTMEs M ( N "' and uncertainties AM N " 
for the (p-p-) Qv decay of 94 ' 96 Zr, 98 ' 100 Mo, 104 Ru, 110 Pd, 
I28,i30rp e an( j 150]^ i so topes. Both bare and quenched values 
of qa are considered. Case I and Case II denote calculations 
with and without SRC1, respectively. 



f3 f3 qa Case I Case II 



emitters 






T7 C U "J 

M N 


AM N 


Mn 


ATT °" 

AM N 


HA r-r 

94 Zr 


1 


254 


126.2146 


44.9489 


152.8378 


27.1912 




1 





142.9381 


49.1752 


172.1620 


29.3965 


96 Zr 


1 


254 


100.5313 


36.8858 


122.5048 


21.9209 




1 





114.4851 


40.3246 


138.6328 


23.5263 


98 Mo 


1 


254 


202.5006 


71.6345 


245.3957 


41.8882 




1 





230.1520 


78.3244 


277.2795 


44.9878 


100 Mo 


1 


254 


206.7533 


73.0792 


250.1870 


43.7119 




1 





235.0606 


79.9885 


282.7964 


47.1334 


104 Ru 


1 


254 


150.5572 


53.9389 


182.7216 


31.9382 




1 





171.8075 


59.0467 


207.1750 


34.3939 


no pd 


1 


254 


231.4743 


82.4924 


280.5688 


49.1588 




1 





263.4339 


90.3033 


317.3947 


53.0150 


128 Te 


1 


254 


126.8285 


46.3381 


153.7370 


29.4676 




1 





143.9772 


50.6942 


173.5263 


31.8554 


130 Te 


1 


254 


136.3856 


46.9164 


164.5378 


27.2226 




1 





154.3797 


51.2511 


185.2849 


29.1907 


150 Nd 


1 


254 


85.5467 


31.4473 


103.4294 


20.9802 




1 





97.3640 


34.5024 


117.0160 


22.8729 



change them by 75%-79%, 58%-62% and 43%-47%, re- 
spectively. Further, it has been noticed through the 
study of radial evolution of NTMEs that the FNS and 



SRC play a more crucial role in the heavy than in the 
light Majorana neutrino exchange mechanism. 

Finally, a statistical analysis has been performed by 
employing the sets of twelve NTMEs M^ v ^ to estimate 
the uncertainties for g^ = 1.254 and = 1.0. It turns 
out that the uncertainties are about 35% for all the con- 
sidered nuclei. Exclusion of Miller- Spenser parametriza- 
tion of Jastrow type of SRC, reduces the maximum un- 
certainties to a value smaller than 20%. The best ex- 
tracted limit on the effective heavy Majorana neutrino 
mass (Mn) from the available limits on experimental 

half-lives Ty 2 using average NTMEs M% v) calculated 
in the PHFB model is > 5.67^94 x 10 7 GeV and 
> 4.061^4 xl ° 7 GeV for 13 ° Te isotope. 
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